Phase Diagram of Electron Systems near the Superconductor-Insulator Transition 
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The zero temperature phase diagram of Cooper pairs exposed to disorder and magnetic field is 
determined theoretically from a variational approach. Four distinct phases are found: a Bose and 
a Fermi insulating, a metallic and a superconducting phase, respectively. The results explain the 
giant negative magneto-resistance found experimentally in In-O, TiN, Bi and high-T c materials. 

PACS numbers: 74.20-z,74.25DW,74.62En 



Introduction. Some alloys such as InO x [HI], TiN [5], 
BiS^Ca^Pri-^Cv^Og+y [6] as well as ultrathin films of 
Bi [7], Be [SHTO] and of high-T c superconductors (SC) 
[6] [11] display (quantum) superconductor to insulator 
phase transition (SIT) (for a recent review see [12]). 
Besides of the very existence of the SIT, the experiments 
demonstrate several unusual phenomena in its vicinity. 
The most general is the maximum of the resistivity 
in external magnetic field in the insulator phase: at 
growing field the resistance first grows, sometimes by 
several orders of magnitude, and then strongly decreases 
[3] [TDJ [13-15 . This surprising decrease is called giant 
negative magnetoresistance (GNM). In [T4UT6] the GNM 
was considered as a signature of localized Cooper pairs 
(CP) surviving in the insulating state and their decay 
in strong enough magnetic field. In the most thoroughly 
studied amorphous alloy InO^ experimenters observed 
a transition from an insulating to a metallic phase at 
increasing magnetic field [3 ], and from superconductor 
to metal if magnetic field and density of carriers grow 
simultaneously [17] • In the same films of InO^, in 
the insulating state at large magnetic field of about 
15T the observed resistivity obeyed Mott's 3d variable 
hopping law R oc exp^o/T) 1 / 4 [14] which shows that the 
Coulomb interaction was screened, and that InO^ films 
of thickness w « 20nra used in [14] must be considered 
as a 3d object. This implies that the essential length 
scales in this alloy were less or of the order of 20nm. 

There are several theoretical approaches to the SIT. 
One of them is based on the BCS theory and accounts 
for the Coulomb interaction enhanced by disorder [18] . 
It shows that in d = 2 even weak Coulomb interaction 
may turn the transition temperature to zero and that 
appearing fermions are localized. Another idea specific 
for granulated superconductors [19] suggests that the 
CP are bound in granules and can tunnel between them. 
Depending on relative strength of Coulomb interaction 
and tunneling amplitude, either the insulating or the 
superconducting state is realized. In [20 the SIT line 
as a function of the CP density n& was calculated for a 
system of charged fermions bound into pairs of a fixed 
size in the field of randomly distributed Coulomb centers. 



In [21] the idea of duality between vortices and CP was 
proposed and some exact results where found. Neither of 
these theories explains the GNM. Numerical calculations 
of the Bogoliubov-de Gennes equations in a random 
medium [22] showed the appearance of superconducting 
islands separated by the insulating see. This islands were 
shown to be gradually suppressed by magnetic field [23] . 

In this article we propose a theory of insulating state 
near the SIT in magnetic field. We assume that the CP 
survive in the insulating state, but are localized by a 
random Gaussian potential. The model also incorporates 
the CP interaction and the magnetic field which destroys 
the pairs when strong enough. Our theory explains 
the GNM and the metal insulator transition (MIT) 
in growing magnetic field and predicts new phases: 
the Bose insulator (BI) and the Fermi insulator (FI) 
[24] . Our physical picture is that in systems displaying 
the SIT the Fermi energy Ep is close to the mobility 
threshold such that kpl « 1. Here I is mean free path in 
the normal state near the SIT. Therefore, the electron 
density fluctuations on large scales like coherence length 
£ and Larkin length Cb(> £ > f° r CP (see below) are 
relatively small. This conclusion is supported by very 
short electron density correlation length found in [22] . 
The disorder does not influence substantially neither the 
number of CP nor their binding energy A. However, the 
density of CP ~ n e A/Ep is much smaller than the 
density n e of the background electrons. Therefore the 
random potential can substantially change the positions 
of centers of mass of CP and localize them. The Coulomb 
interaction at low T seen in experiments is screened on 
distances exceeding the electron distance n e « 7r /ftp, 
which is about lnm [IJG22]). Thus these length scales 
are less than the CP spacing and the film thickness w. 
The Gaussian random short-range potential has a char- 
acteristic localization length C (Larkin length) [25 . The 
random field acting on interacting particles is renormal- 
ized by the self-consistent field and changes with energy. 
The CP and Fermi-excitations appearing at their decay 
are separated by energy gap from the background elec- 
trons. Therefore, the random field acting on the CP and 
excitations is weaker than the field acting on the back- 
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ground electrons. Deeply localized pairs are confined in 
fluctuation potential wells whose size is smaller than C b . 
In such a small well the CP interaction is essential. 
The model and general approach. Let a system of charged 
particles (electrons or CP) in a Gaussian random poten- 
tial U (r) to be exposed to a homogeneous magnetic field 
B. The single-particle Hamiltonian is given by 



U 



2m k 



\2hc) y ' 



+ U k (r), (1) 



where k = b, f for bosons and fermions, respectively, and 
rrif = rri})/ 2 = m and e/ = e b /2 = e. The spin of bosons 
is zero. The random potential has zero average and the 
pair correlator (/7/ c (r)f// c (r / )) = K 2 S(r — r'). The stray 
random potential for electrons is roughly twice smaller 
than that for bosons, i.e. n b « 2k f. Two relevant length 
scales are the magnetic length £ k — \Jhcj (e^B) and the 

Larkin length C k = > /5r(\^7r) d " 2 (ft 2 /(^fe«fc)) 2/(4 " d) J 
where the electron Larkin length C f is 2 4 /( 4 - d ) times 
larger than the CP length C b . At low magnetic fields and 
strong disorder, the CP fill rare deep wells of the random 
potential. We call this state the Bose insulator. When 
the disorder decreases to a critical value, the potential 
wells overlap and transition to superconducting phase 
proceeds. In the insulating phase a strong magnetic field 
destroys the CP by one of two mechanisms: (i) Paramag- 
netic depairing at which Zeeman energy of two electrons 
exceeds the CP binding energy renormalized by disorder 
and (ii) Diamagnetic squeezing at which the size R of 
the optimal potential wells becomes less than the size £ 
of the CP. In a dirty superconductor £ = 0.85^3^/^ 
where £0 = Hvf/(tyA) is the BCS coherence length and 
I is the electron mean free path. Which effect, (i) or 
(ii), dominates depends on d, strength of disorder and 
the electron density. At higher magnetic fields a fraction 
of CPs decays into fermions. As long as the density of 
excited fermions remains small, they are also localized. 
This state we call Fermi insulator. The paramagnetic 
depairing transition from Bose to Fermi insulator (BFT) 
happens when the energy of the CP becomes equal to 
energy of appearing two electrons: 



fj, b -2A = 2(E f +E z ) : 



(2) 



where n b is the chemical potential of the CP at fixed 
density n b and magnetic field B, and 2 A is their binding 
energy; Ef is the disorder energy for an appearing 
fermion and E z = —g e \±BB /2 is its Zeeman energy. 
At further increase of magnetic field the density of 
fermions reaches nf c = Cj d . Then their wave functions 
strongly overlap and their interaction is strong enough 
to overcome the localization provided Cf exceeds the 
Bohr's radius in the media. Derealization of fermions 
turns the FI into a metal, this is the MIT. 
Thin films in parallel field. We first analyze the BFT 
equation for the case of a very thin film of thickness 



iu(<C £,C) in a parallel magnetic field. In this case 
the field does not produce any diamagnetic effect, only 
paramagnetic one. The energy of the first appearing 
fermion Ef can be found from the condition that a 
boson trapped by the same potential well as the fermion 
has energy E b = fi b since levels lower than n b are already 
occupied, whereas E b > \i b does not correspond to 
minimal depairing magnetic field. To calculate Ef we 
apply a variational approach maximizing the probability 
of the potential U(r) at a fixed value of the fermion 
energy Ef\U(v)\ in this potential. The probability P [U] 
for the Gaussian uncorrelated random potential is given 
by P[U]=J DU (x) e"^ f ^ 2 (*)^ Tq find the optimal 
fluctuation (OF) for deeply localized state, it is necessary 
to minimize J U 2 (x) d 2 x at a fixed quantum particle 
energy E U] = f[(H 2 /2m) (V^) 2 + Uip 2 ]dx [28J. The 
minimization over U leads to the relation U = —A?/; 2 , 
where A is a Lagrangian factor. We choose the Gaussian 
fermion trial function ipf(r) = (ctf /7r) 1 ^ 2 e~ afT I 2 . All 
integrals can be calculated explicitly leading to the re- 
sult: A = 2ttH 2 jm and otf = —2mEf/H 2 . Next we solve 
Schrodinger equation for a boson in the potential t/^x) = 
2Uf = —2\i\) 2 p using a similar trial wave function and 
require E b = This procedure gives Ef — (2/9) /i^. 
The chemical potential of the bosons Hb was calculated in 
the works [25 , 27 . For d — 2 it reads n b = — £ b ln(£ b/g^b) 
where £ b = h 2 /(AmC 2 ). The SIT happens at n b g « E b 
where n b w mA/(87rh 2 ). Therefore, only the disorder 
controls closeness to the SIT. From we arrive at the 
critical parallel field £? BFT i| of a 2d film: 
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(3) 



Here B c = 2A/(g e /i#) is the CP breaking field, k = £ b / A 
measures the strength of the disorder relative to the gap; 
its value at the SIT is n c = gn b /A « mg/(87rh 2 ), which 
is the dimensionless (small) CP interaction constant. 

The bosons do not exist at disorder so strong that 
the area of the optimal fluctuation kC 2 / Iii(k/k c ) [25] 
becomes less than 7r£ 2 which happens at n ~ 0.6/ \kpl)- 
For larger k (not shown in Fig. 1) there are no CPs in the 
insulating phase. The squeezing line together with the 
line Q form the BFT line. The fermion density in the 
FI phase is determined by eq. ([3| in which the density 
of bosons must be taken n b = n — rtf/2. At the MIT 
line Ef = —h 2 rif c /(2m) and rif = rif c = Cf~ 2 . Thus 



B 



MIT 



B r 



\ii(k/(k c - ft/32)) 



(4) 



The two lines of the BFT and MIT in the present approx- 
imation cross the SIT line at the same point. The MIT 
line goes to B = 00 at n = 32ft c . Since C b ~ I we rewrite 
Kj ' k c « y/lc/l where (kpl) c ~ 1 at the SIT. A schematic 
phase diagram for the parallel field is shown in Fig. 1. 
Thin films in perpendicular field. In this situation the 
diamagnetic term in the Hamiltonian is nonzero. Still the 
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FIG. 1. Phase diagram for disordered films in parallel fields 
for 7 = 1/2, (kpl) c = 1 and k c — 0.2. The vertical axis 
is the magnetic field in units of the CP breaking field, the 
horizontal axis k F l / (kpl)c- BI, FI, SC and M stand for the 
Bose insulating, the Fermi insulating, the superconducting 
and the metallic phase, respectively. The red curve (online) 
corresponds to the maximum of the resistivity. 



wave functions are isotropic and will be approximated as 



Gaussian ipk = tt 



exp (— akr 2 /2). Otherwise the 



calculation follows the scheme discussed for the parallel 
field case. In this way we find the equation for the energy 
E k = —h 2 ak(l — 3(2/|a/ c ) _2 ))/2m/ c and the density 
of potential wells n w ~ exp[— £ 2 a/c(l — (2a£|) -2 ) 2 ] 
supporting the energy levels lower than a fixed [25] . 
It interpolates between the two limiting cases of energy 
far below and very close to the first Landau level [29] . 
To find the chemical potential of bosons one should 
minimize the total energy /i(a&) = + gnb/n w (Eb) 
including the disorder and interaction contributions over 
a\) (see details in [25 ). For low fields, B < kB c /(1 — 7), 
the results for the case of parallel field apply, here 
7 = 1 — mo/mg e . For larger fields the extension of the 
wave function is essentially determined by the magnetic 
length and the SIT line smoothly crosses over to the up- 
per critical field B c2 = B c d/(8k F l(l - 7)) = c 1 hc/(2e^ 2 ), 
ci w 0.69. At strong magnetic fields, the density n w (E) 
rapidly increases when the energy E^ approaches the first 
Landau level. The value Ef is calculated as before, but 
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FIG. 2. Phase diagram for disordered films in perpendicular 
field using the same parameters and notation as in Fig. 1. 
The BFT line is here due to squeezing. 

with the diamagnetic term. Plugging these values into 



we find the BFT and the MIT lines for the strong 
perpendicular field. Close to the SIT the BFT line is 
given by B^ FT « B^l + ^T -1 - l)«/21n(«/« c )]/7. Far 
from the SIT line the result ([3| applies with the prefactor 
5/18 replaced by 5/9. The squeezing line B sq can be 
more relevant. Its equation reads = £ -2 . This gives 



2B t 



c2 



Cl 



1 — (2cittfcirZln(tt/ ' k c )) 1 



1/2 



(5) 



At 7<2ci/(l + 2ci) and close to SIT (k F l « 1), the 
squeezing transition line passes below the paramagnetic 
one, but deeply in insulator region kpl < 1 they 
interchange. It could explain a controversy between the 
observed strong anisotropy of the resistance in magnetic 
field found in [10] and the absence of anisotropy found 
in [16], both in Be films. The first experiment was done 
very close to the SIT, whereas the second one deeply 
in insulator regime. The phase diagram for this case is 
schematically depicted in Figure 2. 

Three-dimensional case. In d = 3 the directions along 
and perpendicular to B are not equivalent. Therefore, 
the trial wave functions are anisotropic Gaussian, 
(x) = a k p\ ,2 e-^P 2+ ^ It^' 2 . The optimal fluctu- 
ation has ellipsoidal shape with long axis directed along 
magnetic field: /3k — cth — l/(2ttfc^). The squeezing 
transition proceeds when the longitudinal size of the 

— 1/2 

potential well (3 k 1 becomes equal to £. Otherwise 
the calculations are analogous to those in the previous 
sections. The BFT line consists of two parts. The 
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FIG. 3. Phase diagram for bulk samples for the parameters 
of Fig. 1. BI, FI, SC and M stand for the Bose insulating, the 
Fermi insulating, the superconducting phase and the metal- 
lic phase, respectively. This diagram essentially reproduces 
topology of the experimental phase diagram of [12] for In-0 
(right) . The vertical axis of the latter is the magnetic field 
in Tesla, the horizontal axis is again {kpll (hpV) c . 
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squeezing transition obeys equation: 



discussions and B. Roostaei for his help in creating the 
figures. 



B 



3D 



c2 



1 — yj tvkpl ln(ft/ft c 



1/2 



(6) 



where n c ~ mghp/fi 2 is a small number. At decreasing 
disorder, the squeezing curve crosses the paramag- 
netic BFT line whose equation for k ^> k c reads: 
B^ T ^B C [l + 0.2^ ln 2 (ft/ft c )] . For larger densities the 
BFT proceeds due to the Zeeman effect. At n > k c and 
7 < the depairing is caused only by squeezing, eq. 
Equation of the MIT line at k ^> k c is given by B^ T w 

B c 1 + (ft/2) In 2 (ft/ft c ) . For smaller disorder one finds 
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3D 
MIT 



Be 
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(^ln 2 (ft/ft c ) 



1/3' 



Magnetoresistance. In zero magnetic field the wave func- 
tions decay on large scales exponentially, their overlap 
determines the hopping conductivity. For B > the 
decay of the wave function turns into gaussian which 
leads to an increase of the resisteivity [30] . This be- 
havior proceeds until fermions appear in the FI phase 
where the resistance sharply decreases in comparison 
to the BI phase. Indeed, if both resistances have 
variable range hopping nature, they obey the Mott 
formula: R = R exp[- (T /T) 1/(d+1) ]. The differ- 
ence between fermions and bosons is in the value of 
T ~a d / 2 S/n w (/i) - £/ (nC d ). The ratio of the critical 
temperatures reads T //T 6 ~ (n b /nf) (C b / Cf) d+2 . At 
n f ~ n cf = £] 2 one finds T of /T ob ~ (n b /n cb ) (C b /C f ) 2 . 
It is small even in close vicinity of the SIT. In InO^ exper- 
imenters observed the activation behavior of resistance at 
low fields and the VRH behavior at high fields ~ 10— 15T. 
The natural explanation is that the CPs have an intrinsic 
energy gap modified by disorder, whereas there is no gap 
for electrons. 

Conclusions. To conclude, we demonstrated that both 
the BFT and MIT transitions should happen in 2d, 
though the phase diagrams are very different for parallel 
and perpendicular field. We predict a strong anisotropy 
with respect to the direction of magnetic field in 2d close 
to the SIT which must be eventually suppressed by in- 
creasing disorder. 3d systems (films 20nm and thicker) 
are isotropic with respect to the magnetic field direction. 
They always display the BFT transition, but the MIT 
transition does not happen if 7 < 0. Our phase diagram 
in 3d has the same topology as the experimental phase 
diagram established in the review [12] by Gantmakher 
and Dolgopolov. 
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